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sation. Thermal and neutrino-trapping effects are very important ingredients 



of protoneutron stars and the possibility of the delayed collapse during their 
evolution. 
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I. INTRODUCTION 

A protoneutron star (PNS), which is formed after the gravitational collapse of massive 
stars, consists of hot and dense hadronic matter. At birth the matter is lepton rich due to 
the trapping of neutrinos and has entropy per baryon of order 1-3 in units of Boltzmann 
constant. A PNS evolves in a few tens of seconds through the deleptonization and initial 
cooling eras to usual cold neutron stars ||T|. During the evolution, some of the stars may 
become to be gravitationally unstable due to some hadronic phase transitions and collapse 
to the low-mass black holes 0], which we call the delayed collapse. 

To make the low-mass black hole with the gravitational mass of 0(1. 4Mq), first, the 
maximum mass of a PNS must be sufficiently low, and secondly, the mass must exceed the 
maximum-mass during their evolution due to neutrino emission, cooling and/or fall-back of 
matter. The possibility of making the low-mass black holes in supernovae and its relation 
with SN1987A has been discussed by many authors [Q. The observation of neutrinos from 
SN1987A suggested a formation of a PNS in it. However the neutron star has never been 
detected so far. One idea to resolve this mystery is that the neutron star must have collapsed 
further into a low-mass black hole at some later time because of either the accretion of a 
sufficient fall-back mass and/or some change of its equation of state(EOS) 

Thus it is very important to study the change of EOS during the evolution of a PNS to 
discuss the possibility of the delayed collapse: the EOS depends on what kind of hadronic 
phase transitions takes place: quark matter, boson condensed matter, hyperonic matter 
and so on ^. There neutrino-trapping and thermal effects play important roles in the 
rearrangement of constituents of matter and they are also crucial for whether phase transition 
occurs or not. Numerical simulations to describe the evolution of a PNS and the delayed 
collapse have been done [|] 0. In this paper we study kaon condensation in hot and 
neutrino-trapping matter, which may be a key object to the delayed collapse due to the 
large softening of the EOS [0 M. 



Kaon condensation has been studied by many authors mainly at zero temperature [10 , 



since first suggested by Kaplan and Nelson 0. Its implications have been also discussed 
mainly in the context of cold neutron stars; maximum-mass, cooling mechanism, rapidly 
rotating pulsars and so on. Baumgarte et al. |^ performed a numerical simulation of the 
evolution of a PNS by taking into account the change of EOS due to kaon condensation and 
examined how they collapse to black holes. However they used the EOS for kaon condensed 
matter at zero temperature. There has been no study to treat the thermal effect consistently 
in the chiral model. 

Recently we have presented a new framework to treat thermal and quantum fluctuations 
around the condensate on the basis of chiral symmetry ||ll| ||12[ [|13|. Making full use of the 



group structure of chiral manifold, we have derived the thermodynamic potential from the 
nonlinear chiral Lagrangian. The Goldstone mode arises from the V-spin symmetry breaking 
there. We have found the dispersion relation for this mode, and that it plays an important 
role at finite temperature [|14|. Here we study the kaon condensation at finite temperature 



using the chiral model and discuss the properties of the PNS. 

There are many studies about the kaon condensation at finite temperature by using the 
meson-exchange model [|T| ||TB[ . Earlier work suggested that the thermal effect on the EOS is 
relatively small compared with the neutrino-trapping effect ||I||. Pons et al. also discussed the 
EOS for kaon condensed matter and properties of the PNS |]15|. They concluded the thermal 
effect is a key object to induce the delayed collapse. The most important difference between 
the meson-exchange model and the chiral model we use here is about the nonlinearity of the 
kaon field. Since the chiral symmetry dictates the nonlinear Lagrangian for the kaon field, 
the properties of the well-developed kaon condensed phase become very different between 



both models, while those are almost the same as far as the condensate is week [|T| ||T5| . As we 
shall see in this paper the chiral model gives a rather stronger condensation in comparison 
with the meson-exchange model. Thereby the chiral model may give a qualitatively different 
scenario for the delayed collapse of the PNS. 

It is found that a normal PNS without any phase transition cannot collapse under no 
accretion and the delayed collapse is possible by the neutrino-trapping effect in the delep- 



tonization era due to the occurrence of kaon condensation using a most reasonable parameter 
for the KN sigma term. We find that the neutrino-trapping effect gives a main contribution 
to the delayed collapse rather than the thermal effect, within the chiral model. We also 
find, with different parameters, that the thermal effect may be important as well as the 
neutrino-trapping effect. 

In Sect.H we briefly review the essential part of the formulation that enables us to treat 
thermal and quantum fluctuations around the condensate transparently. 

In Sect.HI we discuss some features of EOS in the chiral model and the properties of a 
PNS, especially the possibility of its delayed collapse. 



Summary and concluding remarks are given in Sect. IV. 



II. FORMULATION 

Recently we have presented a framework to take into account the quantum and/or ther- 
mal fluctuations around the condensate, in accordance with chiral symmetry |Tl]| [|T^ ||T3 
The effective partition function Zchirai at temperature T can be written as 

r/3 



^chiral 



N [dU][dB][dB]e-KY> 



^ dTjd^x{CcMral{U,B)+SC{U,B)} 



:i) 



with the imaginary time t = it and P = 1/T. Cchirai{U, B) = Cq{U, B) + CssiU, B) is the 
Kaplan Nelson Lagrangian |^ and is represented with the octet baryon field, B, and the 
chiral field, U = ex'p[2iTa4>a/ f] € SU{3), with SU{3) generators T^ and the Goldstone fields 
(pa- 5i^ is the induced SB term as a result of introduction of chemical potentials, /i^ and /i„, 

- ^tr{i?t[(et[r,^,e] + e[Te™,e1),i?]} 

- :^tr{[Te^, U] [Te^, f/t]} + fXntr{B^B} - /x;^tr{5t[Te„, B]} (2) 

with 17 = ^"^ and Tem = diag{2/3,—l/3,—l/3). Generally it is comphcated to perform 
the path integral by separating the fields (pa into the sum of classical fields {(pa) and 
fluctuation fields (pa in the standard way, because {(pa} reside on the curved manifold 



SU{3) X SU{3)/SU{3). Instead, we introduce the local coordinates around the condensed 
point to eliminate the curvature effect in the neighborhood, which is equivalent to the fol- 
lowing parameterization for f/; 

u = cufa^ = cuy\^ = uuy\), c = ^M^^m/f). (3) 

where (M) represents the condensate, (M) = V+{K+) + V^{K-), with K^ = {(j)i±i(j)^)/y/2 
and 9^ = 2K^ K~ / p , while Uf = exp[2iTa4>a/f] means the fluctuation field. It may be 
interesting to see that this procedure can be also regarded as the separation of the zero- 
mode [|16|. 

Accordingly, defining a new baryon field B' by way of 

B' = u^Bu, (4) 

we can see eventually that 

CcMral{U,B) = Co{Uf,B')+CsB{CUfC,uB'u^), 

6C{U, B) = 6C{CUfC, uB'u^). (5) 

Thus only the SB terms prescribe the KN dynamics in the condensed phase; we can easily 
see that the KN sigma terms, Sxi(^ = ^iP)i stem from Csb: while the Tomozawa- Weinberg 
term from 5C |ll| |13|. The sigma terms can be written by the parameters in the nonlinear 
chiral Lagrangian(ai, 02, 03) and strange quark mass(ms): 

Sicp = -(oi + a2 + 2a3)ms, 

^Kn = -{a2 + 2a3)ms. (6) 

Mass differences of hadrons can determine airug = —67 MeV and 02^5 = 134MeV. There 
remains a large uncertainty associated with the coefficient a^mg. Empirically a^nis lies in the 
range —134 to — SlOMeV, while lattice gauge simulations provide a^nig = —(220 ± 40)MeV 
P^ | |18| ]. Later we mainly use a^rris = — 222MeV as a most reasonable value and add some 



results for the weak limit a^rris = — 134MeV to see the parameter dependence. 



The integration measure is properly approximated as [dUf] ~ Y[a=i [d(pa] in the neighbor- 
hood of the condensate and finally we find 



■'chiral 



\{[d(t)a][dB'][dB']e^V 



a=l 



[^ dTjd'xC:i,JCUf,B') 



(7) 



with the effective chiral Lagrangian, C^jl^j.^i{C,,Uf,B') = Co{Uf,B') + CssiCUfCuB'u^) + 
6C{CUfC,uB'u^). 

We then evaluate the partition function Z^hirai up to the one-loop level under the Hartree 
approximation for kaon-nucleon interactions; we need not care about the nucleon loops in 
the heavy baryon limit (HBL), which we hereafter use in this paper. 

During this course, the poles of the thermal Green function for kaons give the excitation 
spectra of kaonic modes with energies {E±)- The mode with E^ is the Goldstone mode and 
exhibits the Bogoliubov spectrum, which stems from the spontaneous breaking of \^-spin 
symmetry in the condensed phase. It gives a large contribution to the thermodynamic quan- 
tities through the Bose-Einstein distribution function. Under the relevant approximation, 
they are reduced to the simple form 



E±(p) ~ ^p2 + cos(^)m|2 + 52 ± (5 + cos(^)/ii^), (8) 

where m*^ = mj^ — pb/ P(X,KpX + T,Kn{^ — x)) and b = pb{1 + x)/{4:P). We use this form in 



the following calculations. We have discussed its relevance in the previous papers ||T3[ and 
have seen its practical usefulness in giving the thermodynamic quantities. 

Eventually the effective thermodynamic potential fichirai = —T In Z^hirai reads 

^chiral = ^c + ^K + ^A^; (9) 

where Qc is the classical kaon contribution, 

fie = V[-fml{cos{e) - 1) - 1/2 • plf sin2(^)], (10) 

and the thermal kaon contribution, is given as follows; 

Qth^TV /-^ln(l-e-'^^+(P))(l-e-^^-(p)). (11) 
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It is to be noted that the zero-point-energy contribution of kaons is very small [0 [^] and 
we discard it in this paper. 

Qn denotes the nucleon contribution. In order to reproduce the bulk property of nuclear 
matter and get a realistic EOS for kaon condensed matter, we should take into account 
nucleon-nucleon interactions. Following Prakash et al. E^ we introduce the potential energy 
in the form, 

E^/ = E'j^"' + El, (12) 

where £'^™' represents the symmetry energy contribution and E^ the residual potential 
contribution. The symmetry energy is given by 

ET = Vpb(1 - 2xYSP'^\uy, (u = pb/po; Po = 0.16fm-3) (13) 

where the function S^°*(u) reads, 

S^'^'iu) = {So - (2^/3 - l)(3/5)E°)F(n) (14) 

with the constraint F(u = 1) = 1 to reproduce the empirical symmetry energy 5*0 — SOMeV 
at the nuclear density po- Ep is the fermi energy at po- Hereafter we use F{u) = u for 
an example. It is well-known that the nuclear symmetry energy, which stems from the 
kinetic and potential energies, plays an important role for the ground-state properties of the 
condensed phase. Since we have already included the kinetic energy for nucleons, we must 
take into account only the potential energy contribution, which should be beyond chiral 
dynamics. 

On the other hand, as another potential contribution besides the symmetry energy we 
use the following formula which simulates the results given by sophisticated theoretical 



calculations 21 



^X/" - i-^"V« . ^ig^ . 3.,. EC (^J (I - a.cta.|) . (15) 

We choose the parameter set {A,B,B',a,Ci,Ai) for compression modulus to be Kq = 240MeV 



for normal nuclear matter []. Strictly speaking, the effective potential functions Eqs.(|T3]),(|T5D 
should include temperature dependence by way of the Pauli blocking effect on nucleon- 



nucleon scattering in matter. However, as has been shown in ref. [^ temperature- 
dependence should be weak up to several tens of MeV, so that we can use the form E^j^ 
even at finite temperature. 

Including the potential contribution i5^° , we get the single particle energies, 

ep(p) = ^^- i^K, + f^K) (1 - C0S(^)) + -^, 

en(p) = ^- i^Kn + f^K/2) (1 - COs(0)) + -^. (16) 

We have used the non-relativistic approximation for nucleons, where the scalar density is 
not distinguished from the ordinary density pi. 

Then nucleon contributions can be written in the form; 

fi^ = fi^" + Q%°\ (17) 

where ^2^" is the "kinetic energy" contribution and f2^ the "potential energy" contribution. 



Qfeir 

"AT 



2TVy /-^ln(l + e-'^(^'(P)-'^»)), 



<=-PB^. (18) 

OpB 

Here it is useful to introduce the new parameters p^{i = n,p) instead of the chemical 
potentials /ij. 



1^1 = I^P + il^K + Sxp)(l - cos(6')) 



V dp, ■ 



.. , i I^K , ^ \ /-^ „„„//i\\ 1 '^^^N 



^^l = ^^n+[^ + ^Kn) (1 - COs(e)) - -^. (19) 



Then fi^" can be written in the simple form. 



^ Since this parameter set has been also used in Ref. lO] in which V{u) corresponds to our 



expression E^, our results at zero temperature are the same as theirs. 



with the "free" kinetic energy e''(p) = p^ /2M. 
Using the thermodynamic relations, 

Q _ O^ 'chiral ^ _ O^ 'chiral -p — O _l_ T Q -I- X^ C) 

'Jchiral ryp ' t'* a ' -^chiral ^ ^chiral ~r -t ^chiral ~r / ^ fl'i^iy 

we can find charge, entropy and internal energy. The entropy is given by 

Schiral = -f^^chiral + f3V f [-IjfA sin^ {6) + ^xi^ - C0s{6)) 

W) 






-/3\/(1-2x)25^°*(m)P5, 



(20) 



(21) 



(22) 



with the Fermi-Dirac distribution function friE) and the Bose- Einstein distribution function 
fsiE). The kaonic charge is given by. 



Qk = V [i2Kf' sin' {6) + cos(^)n^ + {I + x) sm'{{e)/2)pB 
with the number density of thermal kaons uk, 



Uk 



Nucleon charges are given by 



(2vr)= 



[/B(i?4p))-/B(i^+(P))]- 



Q: = 2V 



d^p 

(2^ 

Then the total energy is given by 

1 
2 
d^p 



/p(e°(p)-^°)l, (^ = n,p). 



E^Mrai = Vf ( --/i^sin^(^) + m^(l - cos(^)) 



+ V 



(27r)= 



[E4p)fB{E_{p)) + i5+(p)/B(i?+(p))] 



+ 2KE/ 



d^p 

pJ (2^ 



eilPJ 



1 (9E 



pot " 



V dpn 



Me°(p)-/i^ 



+ \/(l-2x)^5P°*(M)pB + /ixQi^, 



(23) 



(24) 



(25) 



(26) 



where we add the term with 7-^r , to avoid double-counting. 

V dpn 

The total thermodynamic potential ^totai is given by adding the one for leptons (electrons, 
muons and neutrinos), Vt^ ^totai = ^chirai + ^h 

ni = -2TV J2 I JtS^ [ln(l + e-^(""(P)-'^')) + ln(l + e-^(^^(P)+^'))] , (27) 



'uJ (27r)3 



with £i(p) = e-i{j>) = ^Jmf + p^ {i = e, /x, Ue, v^,). 

The parameters x, {6) and chemical potentials /i;^,/ij(z = n,p) are determined by the 
extremum conditions for given density and temperature. The first condition demands 
dfitotai/dx = 0, which is equivalent with the condition for chemical equilibrium among 
kaons, nucleons: 

l^n- l^p = I^K, (28) 



or 



,0 ,0 ,cPot/.^n n^^ , l-cos(^) 



fi;-fil- ASP°\u){l - 2x) + ^-^^{2{EKn - ^Kp) - fix} + /ii^ = 0, (29) 



by way of Eg . (p!9D . The charge neutrality demands d^totai/dfiR = 0: 

QK + Qe + Q^=Qp{=xVpB), (30) 

where the hadron charges are given in Eqs. (P5|),(P3|) and the lepton numbers are given by 



g, ^ Y,pB = 2^/ ^ [/F(e.(p) - ^.) - /i.(q(p) + p^)] , (31) 

with the lepton- number fractions Yi {i = e, p). Finally the extremum condition with respect 
to {9) gives 

= f smm< - ^f^Kb - /ii cos(^)) + ^^1^, (32) 



which is just the field equation of motion for kaon at momentum zero. Eqs. (29), (30), (32) 
are the basic equations to determine the equilibrated state. 

We can see that once the fluctuation contribution (the second term in Eq.(p2D) is taken 
into account, the relation E_{p = 0) = no longer holds. So the Goldstone nature is violated 

10 



up to one-loop order in the loop expansion. This situation always occurs in the perturbation 
theory [^: the number of loops are always different between the thermodynamic potential 
{Ni) and the self-energy for kaons (A'^; — 1). This shortcoming is never cured unless we 
perform a full order calculation. This matter is beyond the scope of this paper. Fortunately, 
the contribution of fluctuations can be estimated to be small in the temperature-density 



region we are interested in |13 



III. NUMERICAL RESULTS AND DISCUSSIONS 

Using the formulation given in Sect.|I|, we study the properties of kaon condensed matter 
taking into account thermal and neutrino-trapping effects and then discuss the properties 
of a PNS. Finally some implications on its delayed collapse will be given. In this paper, 
we use a^rris = — 222MeV {T,kn = ^{^Kp + '^Kn) = 344 MeV) which is known as a most 
reasonable value (as discussed in Sect.||). In Sect. piID| we will give some results in the case 
of weak coupling limit, a^rris = — 134MeV corresponding to S^'tv = 168 MeV to study the 
dependence of the EOS on the strength of the KN interaction. 

A. Phase Diagram 

First we show the phase diagram in Fig.^ for neutrino-free and trapping matter. In the 
neutrino-trapping case, there exist additional conditions with respect to the conservation 
of lepton numbers [|l7l. Recent calculations of gravitational core-collapse of massive stars 



indicate the electron lepton number Y/g at the onset of neutrino-trapping Yi^, ^ 0.35 p4 
Then we take Y/g = ^e + ^ue = 0-4 and 0.3 and for muon Yi^ = Y^ + Y^, = 0.0 in the 
neutrino-trapping case. On the other hand, in the neutrino-free case, Y^^ = Y^ = 0. 
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FIG. 1. Phase Diagram: solid line and dashed line show the critical lines in the neutrino-trapped 

cases, Yie = 0.4 and 0.3, respectively and dotted line in the neutrino-free case. 

In the phase diagram the temperature-density plane is separated into two regions by the 
critical line; the left-side region means normal phase, and the right-side region the kaon- 
condensed phase. We can see that both of the neutrino-trapping and thermal effects work 
against the phase transition. In the low-temperature case the neutrino-trapping effect is 
dominant, while both effects become comparable at high-temperature. 

The neutrino-trapping effect may be simply understood as follows. Kaon condensation 
occurs when the lowest kaon energy (e;^-) is equal to the kaon chemical potential (/i/^): e^- 
decreases by attractive kaon-nucleon interactions as density increases, while fix = f^n ~ f^p 
increases to match e^- at the critical density. In the neutrino-trapping case, the chemical 
equilibrium due to weak interactions implies that the relation, fix = fJ'u ~ fJ'p = fJ'e — l^v^i 
holds. Compared with the neutrino-free case where iiy^ = 0, fix should be reduced due 
to the finite contribution of /ij,^. Hence higher density is needed to satisfy the condition 
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B. EOS for thermal and isentropic matter 

Here we compare the equations of state in the two cases: the isothermal and isentropic 
cases. The equations of state under the isothermal condition (T = 0, 40 and 80 MeV) 
for kaon condensed matter and normal matter are depicted in Fig.^. The resulting EOS 
for kaon condensed matter shows the existence of the thermodynamically unstable region 
because the phase transition is of first order and large softening of EOS so that we have 
recourse to the Maxwell construction to obtain the realistic EOS in equilibrium^ There 
appears the equal-pressure region as a result of the Maxwell construction. We can see that 
both the thermal and neutrino-trapping effects stiffen the EOS mainly through suppression 
of kaon condensate. Then, in comparison of the EOS with that for normal matter, the 
thermal effect is more profound, in particular around the critical density. 
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^ Strictly speaking, we need to apply the Gibbs conditions because there exist two chemical 
potentials ||2^ |^6|. However, we immediately see that the Gibbs conditions cannot be satisfied in 
the chiral model(Appendix A). So we use here the Maxwell construction for simplicity to get the 
EOS in equilibrium and leave this matter for a future study. 
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FIG. 2. Pressure of isothermal matter (T = 0, 40, 80MeV) for kaon condensed matter[left panel] 
and for normal nuclear matter in beta-equilibrium [right panel]. Solid lines and dashed lines are 
in the neutrino-trapped case (Y^e = 0.4) and the neutrino-free case, respectively. 

The equal-pressure region becomes narrower as temperature increases or more neutri- 
nos are trapped, which implies that the magnitude of the first-order phase transition is 
reduced by both effects of temperature and neutrinos. We can see the remarkable narrow- 
ing, especially in the neutrino-trapping case because trapped neutrinos suppress not only 
the occurrence of kaon condensation but also the growth of the condensate. We shall see 
later that the existence of the equal-pressure region may lead to the gravitationally unstable 
region in the branch of neutron stars. 

The isentropic EOS might be more relevant for the PNS matter. First of all we show 
the isentropic hnes in density-temperature plane in Fig.|| by calculating entropy by the use 



of Eq.(|22D and lepton contribution for given temperature and density. In each isentropic 
line, nucleons dominantly contribute to entropy [|l| P |T^. Under the isentropic condition, 
in which entropy per baryon (S) is taken to be constant (1,2 or 3) over the star interior, 
temperature increases as density becomes high. This means that temperature becomes the 
highest at the center of the star once we construct a PNS with isentropic matter. 
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FIG. 3. Isentropic lines [S = 1,2,3). Solid lines and dashed lines are in the neutrino-trapped 

case (Y/e = 0.4), the neutrino-free case, respectively. 

Then the isentropic EOS can be obtained by connecting the values of pressure at the 
corresponding temperatures for given densities (Fig.H). It is to be noted that the equal- 
pressure region, under the isentropic condition, disappears except the 5* = cases. 
400 
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FIG. 4. Pressure for isentropic matter (S = 0,1,2,3). Solid line and dashed line are in neu- 
trino-trapping case (Y/e = 0.4) and neutrino-free case respectively. 

It is well-known that when the equal-pressure region exists in the EOS, neutron star also 
has a density gap ||2^. In the isentropic case, the equal-pressure region no longer exists and 
the density gap disappears even if we apply the Maxwell construction for the isothermal 
EOS. Instead the mixed phase appears inside the neutron star Q. 

C. Properties of a PNS and the possibility of its delayed collapse 

Using the isentropic EOS, we construct a PNS by solving TOV equation^. In Fig.H we 
show the central density-gravitational mass relation in the neutrino-free or trapping case 
with entropy per baryon 5 = 0, 1, 2 or 3. We can see the thermal effect is large for the stars 
whose central density is around the critical density because the pressure rapidly increases 
in this region as temperature does. On the other hand, mass of the stars with high central 
density seems to be hardly changed by the thermal effect. This is because the thermal effect 
causes not only higher pressure but also larger energy which in turn results in the stronger 
gravity. Sometimes this leads to the lower maximum-mass for the higher entropy. 



^Note that this mixed phase should be distinguished from the one resulted from the Gibbs 
conditions. 



In the low density region Pb < Po, we use the EOS calculated by Lattimer et al [28 
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FIG. 5. Central density versus gravitational mass. Solid lines and dashed lines are in neu- 
trino-trapping (Y/e = 0.4[left panel] or Yi^ = 0.3[riglit panel]) and neutrino-free cases respectively 
at £' = 0,1,2,3. 

We can find that, clearly in the neutrino-free case, once kaon condensation occurs in 
the core of a star, gravitationally unstable region (negative gradient part) appears and 
the neutron star branch is separated into two stable branches: one is for stars with kaon 
condensate in their cores and the other consisting of only normal matter. The maximum 
mass of the stars which include kaon condensate in the core, is hardly changed but decreases 
by the thermal effect because of the strong gravity except the S = 3 case. 

In the S" = 0, 1 and neutrino-trapping cases, we can see that the neutron-star branch 
is also separated by the gravitationally unstable region and the maximum-mass lies on the 
branch with kaon condensate. 

On the other hand, in other cases: the S = 3 and neutrino-free case, the S = 2,3 and 
neutrino-trapping cases, almost all of the stars with kaon condensate are gravitationally 
unstable, and the maximum-mass star still resides on the normal branch. For this reason 
the central density of the maximum-mass star is very different from those in the S* = or 1 
case. In table | we summarize the data for stars with maximum-mass in each case. 
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TABLES 



branch 



central density 



gravitational mass 



total baryon number 



i/-free Yj^^ = 
S = 
S = 1 
S = 2 
S = 3 



K 
K 
K 

N 



13.2/>o 
13.5/90 
12.5po 
4.0po 



1.59M0 
1.59M0 
1.58M0 
1.65M0 



2.24 X lO^'^ 
2.22 X lO^'^ 
2.14 X lO^'^ 
2.11 X lO^'^ 



i^-trapping Yie = 0.4 
5 = K 

5 = 1 K 

5 = 2 N 

5 = 3 N 



10.0/90 

9.6po 
5.4po 
5.3po 



1.64M0 
1.65M0 
1.73M0 
I.86M0 



2.13 X lO^'^ 

2.14 X 10^"^ 
2.23 X lO^'^ 
2.36 X lO^'^ 



i/-trapping Yie = 0.3 
5 = K 

5 = 1 K 

5 = 2 N 

5 = 3 N 



11.3/90 

ll.Opo 
5.1po 
S.lpo 



1.62M0 

1.62M0 
I.66M0 

1.82M0 



2.15 X 10^"^ 

2.14 X lO^'^ 

2.15 X lO^'^ 
2.33 X lO^'^ 



TABLE I. Configuration for stars with maximum-mass: which branch the star resides in, 
normal branch(N) or kaon condensed branch(K), its central density, mass and total baryon number. 
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To discuss the possibility of the delayed collapse of the PNS, the total baryon number 
Nb should be fixed as a conserved quantity during the evolution [^ , under the assumption 
that there is no accretion in the deleptonization and initial cooling eras. 

Discarding gravitationally unstable stars, we show mass of stable stars for given baryon 
numbers in Fig.|^. Each terminal point represents the star with the maximum-mass and 
the maximum total baryon number in each case. If the initial mass is beyond the terminal 
point when the neutron star is born, it should collapse into a black hole (not a delayed 
collapse but the usual formation of a black hole during supernova explosion). We have shown 
the neutrino-trapping and free cases; the former corresponds to the initial stars before the 
deleptonization era, while the latter to the initial cooling era after deleptonization. It is 
interesting to see the difference between the neutrino-free and trapping cases: the curve is 
shortened as entropy increases in the former case, while elongated in the latter case. The 
difference results from which branch the maximum-mass star resides in. These features are 
responsible for the following discussion about the delayed collapse and the maximum-mass 
of cold neutron stars. 
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lines in neutrino-trapping (Yie = 0.4[left panel] or Yje = 0.3[right panel]) and neutrino-free cases, 
respectively at S* = 0, 1, 2, 3. 

The delayed collapse takes place if the initially stable star on a curve finds no correspond- 
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ing stable point on other curves at any stage during its evolution through deleptonization or 
cooling. When we discuss the evolution, several time-scales are important: t^ for the onset 
and growth of kaon condensation, tdei and tcooi for deleptonization and initial cooling eras, 
respectively. We can neglect tx because tx <^ several sec. ~ tdei < 20-30 sec. ^ tcooi |1| 



Now, consider the typical evolution: for example, a PNS has initially YJe = 0.4 and 
S = 2 after supernova explosion and evolves through deleptonization into the S = 2 and 
neutrino- free stage as an intermediate stage. We can see clearly the PNS with I.68M0 < 
Mg < I.73M0 can exist as meta-stable star at the beginning but cannot find any point on 
the curve in the S = 2 and neutrino- free case. Therefore it must collapse to the low-mass 
black hole in the deleptonization era. It is to be noted that because the stars on the curve 
in the S = 2 and Yie = OA case hardly include kaon condensate, their collapse are largely 
due to the appearance of kaon condensate in their cores. 

Furthermore, we can also determine maximum-mass of cold neutron stars by taking into 
account the evolution of the PNS, as pointed out by Takatsuka [^. Usually we assign 
maximum-mass of cold neutron stars as the terminal point on the curve in the S = and 
neutrino- free case in Fig.^ {Mmax = I.SQMq). However, it is wrong when we take into 
account the evolution of neutron stars, especially in the initial cooling era. In order to see 
how to determine the realistic maximum-mass of cold neutron stars, we consider a typical 
evolution with the total baryon number fixed. During the evolution we have already adopted, 
stars with more total baryon number than the terminal point of the curve in the S = 2 and 
neutrino-free case should collapse to black-holes in the deleptonization era and cannot evolve 
stably to the usual cold neutron stars. Therefore only stars with the total baryon number 
Nb < 2.14 X 10^'^ can evolve to cold neutron stars and the corresponding maximum-mass 
is I.54M0 (See Fig.^ and Table |.). Generally speaking, the maximum-mass should be 
determined by taking into account the evolution [^. In our calculation it seems that the 
neutrino-free and hot stage plays an important role to determine the maximum-mass. 

If we consider other scenarios, we can of course determine the maximum-mass of cold 
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neutron stars and stars which should collapse to the low mass blacks hole in each case. 

We can conclude that the delayed collapse possibly takes place in the deleptonization era 
due to the occurrence of kaon condensation and the maximum- mass should be determined 
in the initial cooling era. 

Next, we show that the delayed collapse in the case of the PNS consisting of only normal 
matter (without any phase transition) is impossible. In Fig.^ the gravitational mass for 
gravitationally stable normal neutron stars is plotted as a function of total baryon number. 
Each terminal point represents the maximum-mass star in each configuration. We can 
similarly study the stability of neutron stars in the deleptonization or cooling stage. After 
all every star can find the stable point in each stage, and the delayed collapse of normal 
neutron stars is impossible. 
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FIG. 7. Total baryon number and gravitational mass for normal matter. Solid lines and dashed 

lines are in neutrino-trapping (Y/g = 0.4) and neutrino-free cases, respectively at 5 = 0,1,2,3. 



D. Case of the weak KN sigma term 

The left panel in Fig.|^ shows the isentropic EOS for neutrino-free and trapping (l^e = 0.4) 
matter in the case of "Ekn = 168MeV. The phase transition is of the second order in this 
case. The critical density is moved higher compared to that in the case of T^kn = 344MeV 
(See Fig.^), which leads to the weaker softening of the EOS. There is no thermodynami- 
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cally unstable region in the isothermal EOS and we don't need the Maxwell nor the Gibbs 
constructions. 



400 




3 4 5 6 7 2.35 2.4 2.45 2.5 2.55 2.6 

u= PB /Po Nb (/lO^^ ) 

FIG. 8. The EOS [left panel] and the relation between total baryon number and gravitational 

mass for stable PNS [right panel] in the case of T,kn = 168MeV. Solid lines and dashed lines are 

in neutrino-trapping (Y^e = 0.3) and neutrino-free case at S* = 0, 1, 2, 3. 

With these equations of state, properties of the PNS can be studied as well. The relations 
between the total baryon number and the gravitational mass are shown in the right panel 
in FigJ|. We can also discuss the possibility of the delayed collapse in the case of weak KN 
coupling. The behavior is very different from that in the case of T^kn = 344MeV in the left 
panel in Fig.^. As entropy increases, the terminal point in the neutrino-free case tends to 
move largely to more baryon number and heavier gravitational mass because their central 
densities remain around Gpo, where the thermal effect in the EOS are remarkable in Fig.|^. 
On the other hand, in the neutrino-trapping case, the maximum total baryon number is 
hardly changed in each configuration. 

Then if we assume the same evolution in the above discussion: a newly born PNS 
sitting on the curve in the 5* = 2 and neutrino-trapping case evolves to the corresponding 
one on the curve in the S = 2 and neutrino-free case, delayed collapse can occur in the 
deleptonization era. Moreover, it is also possible in the initial cooling era, different from 
the case of T^kn = 344MeV (See the curves in the neutrino-free case in Fig.^). Therefore 
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we should conclude that the thermal effect may cause the delayed collapse as well as the 
neutrino-trapping effect for the parameter T,kn = 168MeV. 

IV. SUMMARY AND CONCLUDING REMARKS 

In this paper, using our new framework based on chiral symmetry to treat the kaon 
condensation at finite temperature, we have discussed the properties of kaon condensed 
matter and its implication on astrophysics, especially the delayed collapse. 

When we discard any phase transition, the delayed collapse is impossible in deleptoniza- 
tion and initial cooling eras. On the other hand, when kaon condensation is realized, the 
delayed collapse is possible. We found that, in the case of T,kn = 344MeV, which is a 
most reasonable value for the KN sigma term, the neutrino-trapping effect gives rise to the 
delayed collapse under no accretion of matter from the surroundings. This point is a main 
difference to the previous work about kaon condensation and the delayed collapse [0, in 



which the meson-exchange model was used and it was concluded that the thermal effect is 
a key object to the delayed collapse. 

The difference may result from the fact that, in the chiral model, the PNS has very high 
central density, compared with that in the meson-exchange model |l^, due to the stronger 



phase transition. We also studied the case of weak limit of the KN sigma term, S/^jy = 
168MeV, which is a control parameter for the strength of the phase transition in the chiral 
model, to confirm our findings and see to what extent they are modified by the parameter 
dependence. We find that the neutrino-trapping effect can cause the delayed collapse in the 
deleptonization era even in this case. However, we also find that the thermal effect becomes 
important in this limit and the delayed collapse becomes possible by the thermal effect in 
some scenarios. Some people may have impression that the Maxwell construction gives such 
high central density in the chiral model. However, we have rather high central density even 
in the weak coupling limit, as is seen in Sec jlll D] where we no longer need the Maxwell nor 
Gibbs construction. Hence we can say that the nonlinearity of the kaon field, which is an 
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essential feature of the chiral model, gives rise to a stronger phase transition than in the 
meson exchange model. Thus we conclude that the delayed collapse is possible due to the 
occurrence of kaon condensation in the core of the PNS and neutrino-trapping effect gives 
a main contribution in the chiral model. 

In order to study the mechanism of the delayed collapse and mass region of the PNS 
to collapse in more detail, we have to study the dynamical evolution beyond the static 
discussion given here. The accretion of a fall-back mass may be an important ingredient 
P], by which the delayed collapse may occur more easily. Neutrino diffusion is another 
important ingredient in that study: we must treat the change of EOS, neutrino opacity and 
heat capacity in a consistent way which come in the diffusion equation |^ |^ [0. Since 



the neutrino-nucleon scattering or neutrino absorption by nucleons should be enhanced in 
the kaon condensed phase [^ [Q the neutrino-trapping era might last longer if the heat 



capacity is little changed. 

The heat capacity is related to the properties of matter at finite temperature, which 
should be dominated by the nucleon contribution |I|] P [|T^. In the calculations with 



relativistic mean field theory (RMF), the scalar interaction between kaon and nucleon and 
that between nucleons may give contributions through the modification of nucleon mass. 
The effective mass of nucleons affects the entropy or specific heat through the equation for 
nucleon contribution to entropy per baryon (S^), 



S^ ^ xVyMf + (37rp,)2/3 + (1 - x)VVMf + (37rp„)^/3 

7r2T (SttVb)'/^ ^ ^ 



with the low temperature approximation in the relativistic treatment [^ [T^. The nucleon 
effective mass is simply written as M* = M— S/^j(l— cos(6')), {i = n,p) in our case, and there 
is no contribution from the nucleon- nucleon interaction. Then the ratio M* /M remains one 
in the normal phase while it decreases to be 0.6-0.7 in the well-developed kaon condensed 



^The importance of the nonlinearity is also seen in the context of the discussion about the possi- 
bility of the Gibbs construction (see Appendix A). 
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phase. In the reahstic treatments of nuclear matter it should have density dependence even 
in the normal matter, so that it should be a smaller value than we have found; e.g., it 



has been found 0.1-0.4 within the RMF [|T5|. Furthermore we have treated nucleons in a 



non-relativistic way, so that the heat capacity is given only by their kinetic energies without 
any effect by the condensate. Then in the formula of entropy or specific heat the bare mass 
comes in instead of the effective mass and it may lead to the overestimation for specific 
heat. At the same time entropy may be overestimated as well; we may underestimate the 
thermal effect in this paper. In order to improve these points resulting from the nucleon's 
effective mass, we are planning to incorporate the RMF besides the chiral model for the KN 
interaction, as an extension from the previous work at zero temperature p7[| . 

In this paper we studied the kaon condensation and found that the delayed collapse is 
possible due to not only the thermal effect but also the neutrino-trapping effect through 
the suppression of condensate. It is very interesting to study whether the delayed collapse 
occurs in the deleptonization or initial cooling era for other phase transitions, for example, 
hyperonic matter |]^ or quark matter. 
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APPENDIX A 

In this paper we have apphed the Maxwell construction (MC) to get the EOS in equi- 
librium. Strictly speaking, we need to apply the Gibbs conditions (GC) instead of using the 
MC [^ because there exist two chemical potentials: baryon (neutron) and charge (electron) 
chemical potentials. Recently an attempt appeared along this line to derive a realistic EOS 
pi^l . Then EOS should be smoothed due to the appearance of the mixed phase, where normal 
matter and kaon condensed matter coexist. At the beginning of the mixed phase, matter 
is expected to contain droplets of, negatively charged, kaon-condensed matter immersed in, 
positively charged, normal matter. However some works suggest that sometimes the GC 



cannot be satisfied |]15[ @]. Here, we study the behavior of the thermodynamic potential in 
each case of the Maxwell or Gibbs construction, and show that the GC cannot be satisfied 
in the chiral model as already suggested by Pons et al. |jl5[. 

Consider the mixed phase, where the bulk normal matter(A^) and the bulk kaon con- 
densed ma.tteT{K) contact with each other by a sharp boundary^ The MC demands the 
relations 

/^^ = /^^ P'' = P'', (34) 

on the other hand, the GC can be written as 

,,^ ,,K ,,N ,,K r,N r,K /qc\ 

between two phases. Then we can see that the MC corresponds to the Gibbs construction 
in the case of only one chemical potential. First of all we figure out how the MC is applied 
in our study by using thermodynamic potential. 

In Fig.^ thermodynamic potential ^totai/V is shown as a function of density or {6). 
The dotted lines in both panels mean the solutions under the equilibrium conditions; chem- 
ical equilibrium Eq. (l28|) , local charge neutrality Eq. (|30D and equation of motion for kaon 



We discard any effect given by the boundary or the Coulomb energy. 
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at momentuin zero Eg. (|32|) . Then the dotted hnes show the EOS through the relation 
^totai/y = ~P in the left panel, where normal matter with {9) = exists at low den- 
sity pb < Peril = 3.1po and as density increases from pcrit, thermodynamically unstable 
region {dP/dps < 0) appears up to 4.5po and the stable region again appears beyond 
that. On the other hand, the solid lines are calculated with baryon chemical potential fixed 
(/i„ = 1032,1082 and 1132MeV), keeping the conditions for chemical equilibrium Eq.(p8|) 
and local charge neutrality Eq. (|30|) . In this case the cross points of the solid line with 
Pn = 1082MeV and the dotted line in the thermodynamically stable region give the begin- 
ning and ending densities of the equal-pressure region in Fig.H. Hence, if we use the MC 
for the dotted line in the left panel, we get the EOS for the T = and neutrino-free matter 
in Fig.|^. In the right panel, thermodynamic potential in the condensed phase is given as 
a function of the order parameter (6) . The dotted line is given by connecting the extrema 
satisfying dfltotai/d{9) = in solid lines, which corresponds to the dotted line in the left 
panel. It is to be noted that the thermodynamically unstable region in the left panel exactly 
corresponds to the left region to the maximum of the dotted line in the right panel where 
d'^Qtotai/d{0)'^ < 0. The minimum in each solid line means the ground state; in the case 
of low Pn the minimum point stays at (6) = 0, which means the normal phase is favored. 
Afterwards two minima appear at pn = 1082MeV with equal pressure, corresponding the 
coexisting phas^ in the context of the MC. Above this value of baryon chemical potential, 
the kaon condensed phase is always favored. 



^ In this section, we call the mixed phase in the case using the MC as the "coexisting phase" and 
in the case imposing the GC, the "mixed phase" to avoid confusion. 
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FIG. 9. The behavior of thermodynamic potential is shown as a function of density [left panel] 

and {6) [right panel]. /i„ is fixed and chemical equilibrium and local charge neutrality are achieved 

in solid lines. Dotted lines mean the solutions of equilibrium conditions; chemical equilibrium, 

local charge neutrality and equation of motion for kaon. 

Using the MC, typical characteristics are found for the structure of neutron stars as 
follows; 1: Gravitationally unstable region [dMc/dpc < 0) appears. 2: There is no coexisting 
phase in neutron stars^. 3: Equal-pressure region appears in the EOS and it leads to the 
density gap in the structure of a neutron star. 

These features will disappear when we apply the GC, which is a correct treatment when 
there exist more than one chemical potential |]26|. In order to discuss whether the GC 
can be satisfied or not, by the similar manner as in the case of the MC, the behavior of 
thermodynamic potential ^totai/y is shown in Fig.[l^ with two chemical potentials fixed 
at Hn = 1082MeV and /ig = —100, —50,- ■ ■ ISOMeV, keeping chemical equilibrium. If the 
GC can be satisfied, the existence of two minima is expected. However we find only single 
minimum at {$) = 0. Then we can conclude the GC cannot be satisfied in the chiral model 
||r5| . On the other hand, in the meson-exchange model, two minima exist and the GC can be 



applied except the strong phase transition |[T^. The difference between two models results 



As already we noted in Sect. Ill B, this comment is correct only when we construct the neutron 
stars with zero entropy. Finite entropy allows the existence of the coexisting phase. 
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from the nonlinearity of kaon field in the chiral model and we can see that even if we use 
the weak KN sigma term, the GC cannot be satisfied in the chiral modelQ. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

<e> 

FIG. 10. The behavior of thermodynamic potential is shown as a function of {6). Two chemical 
potentials are fixed and chemical equilibrium is achieved in each line. 

It is to be noted that we have discussed the thermodynamic potential in the ideal situation 
where the surface and Coulomb effects are discarded. However, this type of argument may 
miss as essential point for the existence of the mixed phase. When we consider a realistic 
situation for the mixed phase, there is no clear boundary between two phases; the density 
of each constituent and the chiral angle should be smoothly changed there. So the naive 
GC Eg . (^51) cannot be applied to this situation. Instead, we must carefully determine the 
configuration of the mixed phase by allowing the chemical potentials to depend on the spatial 
coordinates and taking into account the screening effect of the electric potential. Then we 
can compare the energy of the mixed phase with the normal phase to see whether the mixed 



^ Using the chiral model and the "potential energy" contribution Eq.([l8|) for A^A^ interaction in 
non-relativistic approximation, the Gibbs conditions cannot be satisfied, even in the case of the 
weak KN sigma term, for example, in the case of Tikn = ISOMeV where appears very narrow 
thermodynamically unstable region. It may also be interesting to study the case of relativistic 
treatment for NN interaction besides the chiral model for KN interaction. 
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phase can exist. Recently there have been appeared some works by taking into account the 
surface and Coulomb effects along this line [Q |3^. In this treatment, we think, there is 
some possibility that the chiral model may give the mixed phase as well. This interesting 
issue is left for a future study. 
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